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INTRODUCTION 

This  document  outlines  a  basic  framework  for  time-domain  solutions  to  Maxwell’s 
equations  in  three-dimensional  geometry.  The  outer  limits  of  the  computational  domain 
consists  of  absorbing  layers  as  described  by  Berenger  (References  1  and  2).  The 
computational  domain  grid  consists  of  square  cells  as  described  by  Yee  (Reference  3). 


MAXWELL’S  EQUATIONS 


Following  Berenger  (References  1  and  2),  we  include  anisotropic  electric 
conductivity  a  and  magnetic  conductivity  <j* .  We  assume  that  the  media  is  non- 
permeable  and  the  permittivity  £  is  non-dispersive.  The  full  set  of  coupled  Maxwell’s 
equations  are  written  as 


c  dt  c 


e  dE„  Atig, 

_ _ Z_ 

c  dt  c 


s  Akc, 

— ^  + ^ 
c  at  c 


e  dEy^  ^  Akg^ 
c  dt  c 


edE 

IX 

c  dt 


+ 


At^ 

c 


(1) 

(2) 

(3) 

(4) 

(5) 
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c  ot  c  '  oy^  '  ^ 


1  dH^  .  47i:&. 


c  dt  c 


1  ,  4;r<T’ 

c  dt 


=-Af£:  +£  ) 
c  a?  c  ^  ax^ 


(6) 

(7) 

(8) 

(9) 


) 


(10) 


(11) 


(12) 


The  purpose  of  anisotropy  in  the  conductivity,  and  the  splitting  of  the  electric  and 
magnetic  fields,  is  to  allow  the  freedom  to  selectively  create  absorption  along  specified 
directions.  Further,  these  conductivities  and  the  field  splitting  are  non-physical  and  are 
introduced  only  to  create  a  perfectly  matched  absorbing  shell  that  surrounds  the 
computational  domain.  The  absorbing  shell  is  designed  to  simulate  an  infinite  sized 
domain  by  virtue  of  the  lack  of  reflected  fields. 

We  now  rewrite  Equations  1  through  12  by  using  identities  of  the  following  form. 
Multiply  Equation  1  by  tx^ATta^tje^  and  we  see  that 

£  d  [■  AitOyt/Sj^  1_  ina^tjs  £  .  ^710^ 

^vj-^  C  dt  C  ^ 
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Comparing  this  with  Equation  1,  we  can  rewrite  Equation  1  (and  similarly,  Equations  2 
through  12)  as 


^|-{exp[4;r(7,r/e]£^}  =  exp[4;r(T//e]£{H„  +  H^] 


c  dt 


dy 


(13) 


f  ^{exp[4;rcT,f/e]^..}  =  -exp[4;rcT,r/e]^{H^,  + 

f  ^{exp[4;r(T,r/e]E^J  =  exp[4;rc7,r/£]^{7/„,  + 

|^{exp[4;rc7//£]£,,}  =  -exp[4;r<T,r/£]|-{/7^  + 

f  ^{exp[4;rcT//£]£„}  =  exp[4;rcT//£]^{i7,,  + 

£  3  3 

f  ^  W^^^/Z^K}  =  -exp[4^<T^V^]^{^^  +  K} 

7I: 

^|:{exp[4;r<7;r]//,J  =  exp[4;r<7;r]^{£,,  +  £,,} 
i^{exp[4;rc7:r]H^}  =  -exp[4;r<T:r]^{E,,  +  £„} 


(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 

(21) 


i^{exp[4;rcT>]/7„}  =  exp[4<r]|-{E„  + 


I 

dx 


(22) 
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^^{exp[4;r(T>]7f^}  =  -exp[4;r(T>]^{£^,  +  (23) 


Equations  13  through  24  represent  our  coupled  set  of  Maxwell’s  equations  to  solve  after 
they  are  written  with  finite  difference  approximations  for  the  space  and  time  derivatives. 
These  equations  are  explicitly  written  with  Cj  ^  0  and  cr*  9^  0  and  this  occurs  only  when 
certain  points  in  the  computational  domain  are  in  an  absorbing  boundary  layer.  In  such  a 
region,  impedance  matching  dictates  that 


In  any  region  of  the  computational  domain  that  is  outside  an  absorption  region, 
(jj  =  <7*  =  0,  and  in  their  continuous  form,  this  means  that  Equations  13  through  24  can 
be  combined  into  the  usual  form  of  Maxwell’s  equations  with 


and  and  E„  +  =  E,  (26) 

and  and  (27) 

SPACE  AND  TIME  DISCRETIZATION 


In  order  to  write  Equations  13  through  24  in  finite  difference  form,  we  let  a  point  in 
space  be  represented  by  (References  3  and  4) 


=  (iAx,7Ay,^Az) 


(28) 


and  a  function  of  space  and  time  by 
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F"  (i,  j,  k)  =  F(iAx,  jAy,  kAz,  nAt) 


(29) 


where  i,  j,  k,  and  n  are  integers.  We  want  to  use  centered  space  and  time  derivatives  to 
achieve  accuracy  0{Ax^),  etc.,  and  O(At^).  Notation  for  space  derivatives  has  the  form 


dF^,k)  _  F^i  + 1/ 2,j,k) - F\i-l/2,j,k)  ^ 


dx 


Ax 


(30) 


and  time  derivatives  are  similarly  written 


dF^Uk)  _  F"^'^2(i,7,fe)-F'-''2M^  +  0(Ar^) 


dt 


At 


(31) 


A  three-dimensional  grid  composed  of  square  cells  will  characterize  the  computational 
domain.  Equations  13  through  24  will  be  discretized  over  a  three-dimensional  Yee  cell 
(Reference  4)  as  shown  in  Figure  1.  For  the  case  where  a  space  grid  point  (i,  j,  k)  lies 
within  the  absorption  shell,  Equations  13  through  24  are  written  differently  for  the  side. 


FIGURE  1 .  Computational  Domain  Unit 
Cell  in  Cartesian  Coordinates. 
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edge,  or  comer  regions  of  the  shell.  In  a  comer  region  of  the  absorption  shell,  we  want 
absorption  to  occur  in  all  three  directions  and  generally,  at  a  given  point  in  a  comer,  all 
three  quWities  are  non  zero  and  different.  In  an  edge  region  of  the  absorption  shell, 

one  of  the  Oj  quantities  will  be  zero  since  absorption  will  occur  in  two  of  three 
directions.  The  remaining  two  nonzero  <7^  quantities  will  generally  be  different  at  a  point 
in  an  edge.  For  either  the  corner  or  edge  region.  Maxwell’s  equations  take  the  general 
form 


FIELD  POINT  IN  CORNER  OR  EDGE  REGION  OF  THE  ABSORPTION 
SHELL 


F;*'(/  +  l/2,7,/:)  =  expl 


' -Anhta^{i  +  \l  2,j,k) 

£(i  +  H2J,k)  J 


EUi  +  \/2,j,k)  + 


cAt 


Ay£(i  +  1/  2,j,k) 


xexp 


f-27cAtcr^(i  +  l/2,j,k)' 
£(i  +  ll2,j,k)  . 


+ 1  /  2,7  + 1  /  2,k)  -  + 1  /  2,7  - 1  /  2,k)  (32) 


+  +  l/2J  +  l/2,k)-  +  \/2,j-l/2,k)} 


£;*'(/  +  !/ 2,7,^)  =  exp 


^  -4nAta^(i  +  l/2,j,k)^ 
£(i  +  l/2,j,k) 


E”  (i +  1/2,  j,k)- 


cAt 


Az£{i  + 1  /  2,j,k) 


xexp 


-27tAtG^  (i  + 1  /  2, 7,  k) 

£(i+T72jJ)  ^ 


+  l/2,j,k  +  l/2)-  + 1  /  2,j,k  -1/2)  (33) 


+  + 1  /  2,i,k  + 1  /  2)  -  +  \l2,j,k-\l2)} 
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+ 1  /  2,/:)  =  expf  +  lEliiJ  +  l/2,k)  +  - — 


£{i,j  +  \/2,k) 


Az£(i,j  +  11 2,  k) 


-2nAta^ {i,i +  l  1 2,k)\ ^ .  . 


X  exp  \{K  +  l/2,k  +  U2)-  +  l/2,k-l/2)  (34) 

V  "f"  1  /  a/c)  j 


+  +  l/2,k  +  l/2)-  +  \l2,k-\l2)} 


+  l/2,k)  =  expf  J  + 1  >  2-,k)  U  ^ ^  j  /  2,jt) - - 

^  \  e(i,j  +  l/2,k)  j  Ax£(/,j  +  1/2,^) 


—2icAt(j^{i,j  +  \  1 2,k) 
E{i,j +  \  1 2,k) 


\i  + 1  /  2,7  + 1  /  2,k)  -  - 1  /  2,;  + 1  /  2,k)  (35) 


+  + 1  /  2,7  + 1  /  2,k)  -  - 1  /  2,7  + 1  /  2,k)] 


E"*\i, i,k  +  \  1 2)  =  expf  ’f  ^  (i, 7, /:  + 1  /  2)  +  - — 


£{i,j,k  +  \l  2) 


Ajc£(/,7,/:  +  1/2) 


X  expf  2^g.(r,7^^  1 1  /  2)  + 1  /  2,7,^:  + 1  /  2)  -  - 1  /  2,7,*  + 1  /  2)  (36) 

V  £(j,7,/:  +  1/2)  /  ■ 


+  + 1  /  2,j,k  + 1  /  2)  -  - 1  /  2,7,A:  +  1/2)} 
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,  f  -47cAt<y.. (i,j,k  +  l/  2)'\ 

E”:\i,j,k  +  1/2)  =  exp( - E"^ii,j,k  + 1  / 2) ■ 


€ii,j,k +  1/2)  ) 


edit 


Aye{i,i,k +  1/2) 


xexp 


^-27iAt(jy{i,i,k  +  l/2)\^^ 


\{H':^'\i,j  +  l/2,k  +  l/2)-H':^'\i,j-l/2,k  +  l/2)  (37) 
e{i,i,k +1/2)  )'■  ■ 


+  +  l/2,k  +  l/2)-  -l/2,k  +  l/2)} 


+  l/2,k  +  l/2)  =  txp(-47tAt<J*ii,j  +  l  /2,k  +  l/  2))Hi:''\i,j  +  l/2,k  +  l/2) 

-  ^exp{-27tAt<7l(i,j  +  l/2,k  +  l/  2)){EliiJ  +  l,k  +  \/2)-  E"Ji,j,k  +  1/2)  (38) 

+ E;(i,j +i,k+i/2)-  E;(i,j,k +1/2)} 

+  l/2,k  +  l/2)  =  exp{-4nAtal(i,j  +  l/2,k  +  l/  2))H'';''\i,j  +  l/Xk  +  l/2) 

+  — exp(-2;iAr<7;(j,;  +  l/2,k  +  l/  2)){E’'^{i,j  +  l/2,k  +  l)-E1{i,j  +  l/  2,k)  (39) 

Az 

+  E^iUj  +  ll2,k  +  l)-  E;^(i,j  + 1  /  2,k)] 

+ 1  /  2,j,k  + 1/2)  =  cxp{-47rAtal(i  +  l/2,j,k  +  l/  2))H;-^'\i  +  l/2,j,k  +  l/2) 

-  ^txp{-2nAtal(i  +  l/2,j,k  +  l/  2)){£;(j  +  l/2,j,k  +  l)-  E^d  + 1  /  2,j,k)  (40) 

+  E"Ji  + 1  /  2,j,k  + 1)  -  E"Ji  + 1  /  2,j,k)} 
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+ 1  /  2,i,k  + 1  /  2)  =  exp(-4Mf(T;(/  + 1  /  + 1  /  2))H;;''\i  +  \  /  2,j,k  +  l  /  2) 

+  ^exp{-2nAtal(i  + 1  /  2,j,k  + 1  /  2)){£;(/  + 1,;,/:  + 1  /  2)  -  El(i,j,k  +  1/2)  (41) 

+  E;{i  +  lJ,k  +  U2)-  E;(i,j,  A:  + 1  /  2)} 


+ 1  /  2,;  + 1  /  2,ik)  =  t\p(-A7cAt(Tl(i  +  l/2,j  +  \/  2,k))H^j''\i  + 1  /  2,;  + 1  /  2,k) 
-^txp{-2nAialii  + 1  /  2,;  + 1  /  2, /:)){£;</  +  IJ  +  \/2,k)-  E"Ji,j  +  l/2,k)  (42) 

+  £;(/  + 1,7  + 1  /  2,A:)  -  E^iiJ  + 1  /  2,  A:)} 


+ 1  /  2,7  + 1  /  2,  A:)  =  exp(-4Mr(7;(j  + 1  /  2,7  + 1  /  2,A:))//;:'''(/  + 1  /  2,7  + 1  /  2,k) 

+  ^exp{-27tAtCT;(i  + 1  /  2,7  + 1  /  2,A:)){£;(i  + 1  /  2,7  +  1,A:)  -  E^d  + 1  /  2,7,A:)  (43) 

+  £:;(/  + 1  /  2,7  +  1,A:)  -  E"Ji  + 1  /  2,7,A:)} 


In  the  event  a  point  lies  in  a  side  region  of  the  absorption  shell,  then  two  of  the  (Jj 
quantities  vanish  and  the  remaining  Gj  quantity  creates  absorption  in  a  direction 

perpendicular  to  the  side.  In  this  case,  by  using  Equations  26  and  27,  Equations  32 
through  43  simplify  from  12  unknowns  to  10  unknowns. 
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FIELD  POINT  IN  ABSORPTION  LAYER  THAT  IS  PERPENDICULAR  TO 
X  DIRECTION:  (7,  =  cr^  =  0 

£\J,  +  1  /  Z,J,K) 

+  l/2J-l/2,k)  +  +  U2,j  +  l/  2,k)  -  + 1  /  2,;  - 1  /  2,k)]l^y 

(44) 

+ 1  /  2,  j,jk  + 1  /  2)  -  + 1  /  2,i,k  - 1  /  2)  +  + 1  /  2,j,k  +  1/2)- 

H;:^'\i  +  \l2,i,k-\l2)}l^z] 


+i/2,k)  =  E;^(i,j +i/2,k)+ 


cAt 

£(/>/  +  1/2,/:) 


+ 1  /  2,)k  + 1  /  2) -  +  l/2,k-l/  2)]/Az 


(45) 


E";\i,j  +  l/2,k)  =  exp 


^  -AnAto^  {i,j  +  \l2,k)^ 
^  e{i,j  +  \l2,k)  ) 


(/,7  +  l/2,/:)  + 


cAt 

e(/,j  +  1  /  2,k) 


(  27rAt<j^(iJ  +  l/2,k)  + 1  /  2, 7  + 1  /  2, it)  -  //"*'''(/  - 1  /  2, ;  + 1  /  2, /:)  (46) 

1,  e(j,;  +  l/2,/:)  “ 

+  l/2,j  +  l/2,k)-  - 1  /  2,;  + 1  /  2,/:)}/Ax 
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E"*\iJ,k  +  l/2)  =  exp\ 


f  -4;rAtcrJiJ,k  +  1  /  2)) 


V  e(iJ,k  +  l/2)  j 


E^(U,k +  1/2)  +  — 


cAi 


e(i,J,k  +  l  /  2) 


xexp 


-2jtAt<7ji,j,k  +  l/2)]f 


+ 1  /  2)  -  - 1  /  2,7,*  +  1/2)  (47) 

£(i,7,*  +  l/2)  y*- 


+  + 1  /  2,7,*  + 1  /  2)  -  - 1  /  2,7,*  + 1  /  2)}/Ax 


rAt 

E:;\i,j,k + 1  / 2) = £;(i,7,* + 1  / 2) - 
^  e(i,7,*  +  l/2) 


X  + 1  /  2,*  + 1  /  2)  -  - 1  /  2,*  + 1  /  2)}/Ay 


H"  (i,j  + 1  /  2,*  + 1  /  2)  =  + 1  /  2,*  + 1  /  2) 


(48) 


•  cA?[{£;(i,7  + 1,*  + 1  /  2)  -  £;(i,7,*  + 1  /  2)  +  ^;(/,7  + 1,*  + 1  /  2) 

-  £;(/.7,A:  + 1  /  2)}/A>;  +  {£;(/,7  + 1  /  2,*  + 1)  -  E;ii,j  +  1/2,*) 

+  E^iiJ  + 1  /  2,*  + 1)  -  £;(i,7  + 1  /  2,*)}/Az] 


(49) 


//"*‘'^(/  + 1  /  2,7,*  + 1  /  2)  =  + 1  /  2,7,*  +  1/2) 


-  cAt{E"ii  + 1  /  2,7,*  + 1)  -  £:;(/  + 1  /  2,7,*)}/Az 


(50) 
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+ 1  /  2,j,k  + 1  /  2)  =  exp(-4M/cT;(i  +  \l2,j,k  +  \l  2))H;:''\i  + 1  /  2,j,k  +  1/2) 

+  cA?exp(-2;rA/<T*(/  +  1  /  2,j,k-¥\l  2)){£’^(j  +  l,j,^  +  1  /  2)  -  £^(/,j,/:  +  l  /  2)  (51) 

+£;(/  + 1,;,/:  + 1  /  2)  -  Elii,i,k^  1  /  2)}/Ax 

+  l/2,j  +  l/2,k)  =  exp(-4;tArcr*(«  + 1  /  2,;  + 1  /  2,k))H"~''^  {i  + 1  /  2,;  + 1  /  2,A:) 

-  cAr  exp(-2;rAr(7*  (i  + 1  /  2, ;  + 1  /  2,  /:))•[£"(/  + 1, ;  + 1  /  2,  A:)  -  (/, ;  + 1  /  2,  ^)  (52) 

+£;(/  +  \,j  + 1  /  2,/:)  -  £:;(i,;  + 1  /  2,/:)}/Ax 

//;;‘''(t + 1  /  2,  j  + 1  /  2,/:)  =  //;■'''(/ + 1  /  2,7 + 1  /  2,k) 

(53) 

+  c^t{El{i  + 1  /  2,;  + 1,/:)  -  Ey  + 1  /  2,;,/:)}/Ax 


FIELD  POINTS  IN  ABSORPTION  LAYER  THAT  IS  PERPENDICULAR  TO 
y  DIRECTION:  =  (t,  =  0 


NAWCWD  TP  8447 


+ 1  /  2,j,k)  =  E”Ji  + 1  /  2,j,k)  — - 


e(i  +  \/2,j,k) 


X  + 1  /  2,j,k  + 1  /  2)  -  + 1  /  2,;,/:  - 1  /  2)}/Az 


+ 1  /  2,/:)  =  E"(i,j  + 1  /  2,^)  + 


£(i,j  +  l/2,k) 


+  l/2,k  +  l/2)-  +  \l2,k-\l2)  +  +  l/2,k  +  l/2) 

( 

+  l/2,k-l/  2)]/Az  +  + 1  /  2,7  + 1  /  2,*) -  1 2,] +\  1 2,k) 

+  + 1  /  2,7  + 1  /  2,fc)  -  - 1  /  2,7  + 1  /  2,^)}/Ax] 


£;^'(/,7,A:  + 1  /  2)  =  El{i,j,k  + 1  /  2)  +  -— 


£({,7,/:  +  1/2) 


X  + 1  /  2,7,  A:  + 1  /  2)  -  - 1  /  2,7,  A  + 1  /  2)}/Ax 


»*i  ,  f -4;zA?<Tji,7,A:  +  l/2))  „  cA? 

r;'(i,7, it  + 1  /  2)  =  exp  - ^  \El{i,j,k  + 1  /  2) - 

I  £(j,7,A:  +  l/2)  )  ^  £{i,j,k  +  \l2) 


( -2n^tcAi,j,k  +  \l  2)\ 

X exp  -  7^  - P;  0'.;  + 1  / 2,A  + 1  / 2) - A/;-^''^(/,7 - 1  / 2,/:  + 1  / 2)(58) 


+ 1  /  2,A:  + 1  /  2)  -  iy;;"''^(i,7  - 1  /  2,  A  + 1  /  2)}/A3; 
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+ 1  /  2,k  + 1  /  2)  =  e\p{-47rAia*XiJ  +  l/2,k  +  l/  j  +  1/ 2,^  +  1/ 2) 

-  cAtexp(-27cAtal(i,j  +  l/2,k  +  l/  2)){El(i,j  +  l,k  +  l/2)-  E"Ji,j,k  +  1/2)  (59) 

+  E;(i,j  +  l,k  +  l/2)-  E;Xi,j,k  + 1  /  2)}|^y 


+ 1  /  2,ife  + 1  /  2)  =  +  l/2,k  +  \/2) 

+  cAt{E;(i,j  +  l/2,k  +  l)-  E;{iJ  + 1  /  2,k)}/Az 


(60) 


+ 1  /  2,7,ik  + 1  /  2)  =  +  l/2,j,k  +  l/2) 

-  cAr[{£:;,(/  +  l/2,j,k  +  l)-  E^Xi  + 1  /  2,j,k)  +  £;(/  +  !/  2,7,/:  + 1) 
-£;(/  + 1  /  2.7,/:)}/At  -  {£■;(/  + 1,7,/:  + 1  /  2)  -  El^{i,i,k  + 1  /  2) 
-£;(/  + 1,7^^  + 1  /  2)  +  £;(i,7\^  + 1  /  2)}/Ax] 


(61) 


^«*i/2( .  + 1  /  2,7  + 1  /  2,  Jfc)  =  //;■■''(/  + 1  /  2,7  + 1  /  2,/:) 
-  c^t{E;  a +ij+i/2,k)-E;  (/,  7+1/2,  k)}/Ax 


(62) 
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+ 1  /  2,7  + 1  /  2,k)  =  exp(-4;zAffT;(z  + 1  /  2,7  + 1  /  + 1  /  2,7  + 1  /  2,k) 

+  cMtX'p[-2JtMc]{i  + 1  /  2,7  + 1  /  2,ifc)){£:;(i  + 1  /  2,7  + 1,^)  -  £;(z  + 1  /  2,7,/:)  (63) 

+  £;(z  + 1  /  2,7  + 1,/:)  -£;(/  + 1  /  2,7,^)}/Aj 


FIELD  POINTS  IN  ABSORPTION  LAYER  THAT  IS  PERPENDICULAR  TO 
Z  DIRECTION:  <t^  =  =  0 


+ 1  /  2,i,k)  =  Eld  + 1  /  2,7,*)  + 


cAt 

£(Z +1/2,7,*) 


X  {//;"''^(z  + 1  /  2,7  + 1  /  2,*)  -  "^(z  + 1  /  2,7  - 1  /  2,k)}/Ay 


(64) 


£”'(z  +  l/2,7,*)  =  exp 


(  ~4Mrc7^(/  + 1/2,7,/:) 
e(z +  1/2,7,*) 


r  (z  + 1/2,7,*)- 


cAf 


f(z  +  l/2,7,*) 


xexp 


-2;rAf<7,(z  +  l/2,7,*)V„ 


wo  7;  {*/r^(^  +  l/2,7.*  +  l/2)-Hr'^(z  +  l/2,7,*-l/2)  (65) 

£(z +  1/2,7,*)  J 


+ 1  /  2,7,*  + 1  /  2)  -  //;*'''(/  + 1  /  2,7,*  - 1  /  2)}/Az 


V  £(z,7+1/2,*)  y 


r(z,7  +  l/2,*)  +  -— 


cAz 


£(z,7  + 1  /  2,*) 


xexp 


-2Mzt7^(z,7  +  l/2,*)V„, 


'  + 1  / 2,*  + 1  / 2) - //;*''^(z,7  + 1  / 2,* - 1  / 2)  (66) 

^(^>7  +  1/2,/:)  y 


+ 1  /  2,*  + 1  /  2)  -  + 1  /  2,*  - 1  /  2)}/Az 
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X  {H";''\i  + 1  /  2,7  + 1  /  2,it)  -  - 1  /  2,7  + 1  /  2,/:)}/Ax 


(67) 


E'':\iJ,k  + 1  /  2)  =  £;(/,7,A:  + 1  /  2)  +  — 


cA? 


e(/,7,A:  +  l/2) 


X  [{^r''(*  + 1 ' 2.7.^  + 1  /  2)  -  - 1  /  2,7,/:  + 1  /  2) 


+ 1  /  2,7,/:  + 1  /  2)  -  - 1  /  2J,k  + 1  /  2)}/Ax 


(68) 


■  {7/;:'''(/,;  + 1  /  2,;t  + 1  /  2)  -  (/,7  - 1  /  2,/:  + 1  /  2) 


+  + 1  /  2,/:  + 1  /  2)  -  /7;""^(/,7  - 1  /  2,/:  + 1  /  2)}/Ay] 


+  l/2,k  +  l/2)  =  + 1  /  2,/:  + 1  /  2) 

-  cAr{^;(/,7  + 1,/:  + 1  /  2)  -  £;(i,7,/:  + 1  /  2)}/Ay 


(69) 


+ 1  /  2,/:  + 1  /  2)  =  exp(-4;rA?(T*(t,7  + 1  /  2,/:  + 1  /  2))//”"‘'^(i,7  + 1  /  2,/:  + 1  /  2) 

+  cA/exp(-2M?o'*(/,7  + 1  /  2,/:  + 1  /  2)){£”j(i,7  + 1  /  2,k  + 1)  -  + 1  /  2,^)  (70) 

+  E^iiJ  +  U2,k  +  \)-  E-JiJ  + 1  /  2,k)}/Az 


16 


NAWCWD  TP  8447 


+ 1  /  l,j,k  + 1  /  2)  =  exp(-4;rAfcr;(/  + 1  /  2,7,it  + 1  /  2))H"-''\i  + 1  /  2,;,  A:  + 1  /  2) 

-  cA/exp(-2MfcT;(/  + 1  /  2,7, ^  + 1  /  2)){£;(/  + 1  /  2,;,^  + 1)  -  £;(i  + 1  /  2,;,^)  (71) 

+  E-^ii  + 1  /  2,7,A:  + 1)  -  £”  (/  + 1  /  2,7,fc)}/Az 


+ 1  /  2,;,*  + 1  /  2)  =  + 1  /  2,;,/:  + 1  /  2) 


+  c/St {eI (i  +  \,j,k +  1/2)-  El (i,  j,k  +  \l  2)] /Ax 


(72) 


//f '"(/  + 1  /  2,7  + 1  /  2,A:)  =  //;■’''(/  + 1  /  2,7  + 1  /  2,*) 

-  cAr[{£;(l  + 1,7  + 1  /  2,^)  -  r  (1,7  + 1  /  2,A:) 

+  £”  (i  + 1,7  + 1  /  2,^:)  -  ^;(j,7  + 1  /  2,^)}/Ax  (73) 

-  {£;(/  + 1  /  2,7  +  U)  -  Eld  + 1  /  2,7,/:) 

+  E"Ji  + 1  /  2,7  + 1,/:)  -  £;(/  + 1  /  2,7,/:)}/A3;] 


However,  since  most  of  the  computational  effort  will  be  spent  outside  an  absorption 
region,  the  Berenger  anisotropic  absorbing  conductivities  (7j  and  a*  (j  =  x,y  ,  z)  will  be 

set  to  zero  most  of  the  time.  In  this  case.  Equations  32  through  43  will  simplify.  If  we 
are  outside  of  a  Berenger  absorbing  region,  but  there  are  physical  regions  consisting  of 
ordinary  conducting  media  in  the  computational  domain,  rather  than  set  (jj  and  a*  =  0, 

we  can  let  Oj  =  a  and  cr*  =  0  where  <T  is  a  spatially  dependent  scalar  conductivity.  In 

the  event  we  are  outside  of  a  Berenger  absorbing  region  and  there  are  no  conducting 
media  (c  =  0),  then  the  spatial  dependence  of  e  describes  the  dielectric  permittivity  which 
characterizes  the  structure  being  studied  (waveguide  or  whatever).  For  the  case  where 
(jj  =  <7,  and  <7*  =  0,  Equations  32  through  43  simplify  to 
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FIELD  POINT  OUTSIDE  ABSORPTION  SHELL:  <7^  =  cr,  =  (T,  =  <7  AND 
<7:=c7;=(t:=o 


E"*‘(i  +  l/2J,k)  =  exp 


-4nA,a(:  + 1  /  2.M)  1  /  2,j.k)  +  — ^ 


e(i  +  l/2J,k)  ) 


e{i  +  \l2,j,k) 


X  exp 


-2nbda{i  + 1  /  2,;,/:) 
£0-+i/2,j,*)  ; 


Hr'\i  + 1  /  2,;  + 1  /  2,ifc)  -  r  (/  + 1  /  2,;  - 1  /  2,  A:) 

Ay 


(74) 


+ 1  /  + 1  /  2)  -  + 1  /  2,;,/:  - 1  /  2) 
hz 


F;"'(fJ  +  l/2,*)  =  exp 


^-4jtAta(i,j  +  \/2,k)^ 

£(i,j+i/2,k)  ; 


£;(z,;  +  l/2,A:)  +  — 


cAr 


e(i,j  +  l/2,k) 


(- 


X  exp 


■2nAta{i,j  +  1!  2,  k) 


y  £{Uj^\l2,k)  J 


+  \l2,k  +  \l2)-  +  \l2,k-l/2) 

Az 


(75) 


+ 1  /  2,7  + 1  /  2,k)  -  - 1  /  2,;  + 1  /  2,^)' 

Ajc  ^ 


£r'a,y,A:  +  l/2)  =  exp 


~4;rArcj(f,  7,  A:  + 1  /  2) 
£(z,7,/:  +  1/2) 


£;(z,;,A:  +  l/2)  +  — 


cAt 


e{i,i,k  +  \l2) 


f . 


X  exp 


-27cAt(j(iJ,k  + 1  /  2)  "if  + 1  /  2.;,^  + 1  /  2)  -  - 1  /  2,7./:  + 1  /  2) 


£(,ij,k  +  ll2)  J 


Ax 


(76) 


+ 1  / 2,it  + 1  /  2) - - 1  /  2,k  +  l/2)  [ 


Ay 
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i/;  (j,  J  + 1  /  2,ik  + 1  /  2)  =  +  l/2,k  +  \/ 2) 


E"(i,j  +  l,k  +  l/2)-  E”(i,j,k  +  1/2)  _  +  \/2,k  +  l)-  E^jiJ  + 1  /  2,^)] 

[Ay  Az 


+  l/2J,k  +  l/2)  =  + 1  /  2J,k  +  1/2) 


f  £;(/  +  l/2,j,k  +  l)-  Elii  + 1  /  2,;,/:)  £;0-  +  \,j,k  + 1  /  2)  -  +  1/2)^ 

Az  Ax 


(77) 


(78) 


''(/  + 1  /  2,  j  + 1  /  2,it)  =  + 1  /  2,;  + 1  /  2, it) 


'£"(/  + 1,7  + 1  /  2,/:)  -  £:;(/,7  + 1  /  2,/:)  £"(/  + 1  /  2,7  +  l,it)  -  EUi  + 1  /  2, 7, it)' 

~cArs  - - ^ 

Ax  Ay 


(79) 
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